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We study the strangeness electromagnetic form factors of the nucleon from the Nf — 2 -\- 1 
clover fermion lattice QCD calculation. The disconnected insertions are evaluated using the 
Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expan- 
sion. In addition to increasing the number of Z(4) noises, we find that increasing the num- 
ber of nucleon sources for each configuration improves the signal significantly. We obtain 
G\j{Q) = —0.017(25) (07), where the first error is statistical, and the second is the uncertainties 
in and chiral extrapolations. This is consistent with experimental values, and has an order of 
magnitude smaller error. We also study the strangeness second moment of the partion distribution 
function of the nucleon, {x^)s-s- 
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1. Introduction 

In the pursuit of the full understanding of QCD, it has been essential to study the structure 
of the nucleon. The strangeness content of the nucleon particularly attracts a great interest lately. 
As the lightest non-valence quark structure, it is an ideal probe for the virtual sea quarks in the 
nucleon. Recently, intensive experiments have been carried out for the electromagnetic form fac- 
tors by SAMPLE, A4, HAPPEX, GO, through parity-violating electron scattering (PVES). The 
global analyses [|l||, i have produced, e.g., G^e^) = -0.008(16) and G^(e^) = 0.29(21) at 
q2 = 0. IGeV^ |@], but substantial errors still exist so that the results are consistent with zero. Mak- 
ing tighter constraints on these form factors from the theoretical side is one of the challenges in 
QCD calculation. Moreover, such constraints, together with experimental inputs, can lead to more 
precise determinations of various interesting quantities, such as the axial form factor G\ [^], and 
the electroweak radiative corrections including the nucleon anapole moment, Ga Unfortu- 
nately, the theoretical status of strangeness form factors remains quite uncertain. For instance, the 
values for the magnetic moment G^(0) from different model calculations vary widely, from —0.5 
to +0.1. One may put an expectation on lattice simulation as the most desirable first-principle 
calculation. The evaluation of strangeness matrix elements, however, has been a serious challenge 
to the lattice QCD, since it requires the disconnected insertion (DI), for which the straightforward 
calculation requires all-to-all propagators and is prohibitively expensive. Consequently, there have 
been only few DI calculations, where the all-to-all propagators are stochastically estimated, and 
quenched approximation is used with Wilson fermion [||, ^ ^. There are also several indirect 
estimates using quenched [||, ||] or unquenched [10| lattice data for the connected insertion (CI) 



part, and the experimental magnetic moments (or electric charge radii) for octet baryons as inputs. 
In this proceeding, we present the first full QCD lattice simulation of the direct DI calculation for 



strangeness form factors. For details of this study, see Ref . [11] 



We also study the strangeness contribution to the lowest moments of the parton distribution 



function of the nucleon, (x), (x ) [12]. In particular, the strangeness second moment {x )s^s 



J dxx^{s{x) — ^(x)) is important quantity in relation to the NuTeV anomaly. The NuTeV experi- 
ments [ p3| ] measured the Weinberg angle, which is 3 a away from the world average value. While 
this may indicate the existence of New Physics, it is necessary to make exhaustive investigations on 
hidden systematic uncertainties before making such a conclusion. In fact, the asymmetry between 
strange and anti-strange parton distribution is one of the probable candidates to explain the NuTeV 
anomaly without New Physics. We note that lattice QCD can make a great contribution to this 
issue, by measuring the strangeness asymmetry in the second moment, which provides essential 
constraint to the asymmetry in the strangeness parton distribution function. In this proceeding, we 
present a preliminary result of (x^)^-, as the first full QCD calculation for this quantity. 

2. Formalism and lattice calculation parameters 

We employ A^y = 2 + 1 dynamical configurations with nonperturbatively i&{a) improved clover 



fermion and RG-improved gauge action generated by CP-PACS/JLQCD Collaborations [14]. We 
use j3 = 1.83 and Csw = 1.7610 configurations with the lattice size of X r = 16^ X 32. The lattice 
spacing was determined as a^^ = 1.625GeV, using A'-input or 0-input [11^]. For the hopping 
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parameters of u,d quarks (Kud) and s quark (Kg), we use K"„^ = 0.13825, 0.13800, and 0.13760, 
which correspond to rriji = 0.60, 0.70, and 0.84 GeV, respectively, and = 0.13760 is fixed. 
We perform the calculation only at the dynamical quark mass points, where 800 configurations 
are used for = 0.13760, and 810 configurations for = 0.13800, 0.13825. The periodic 
boundary condition is imposed in all space-time directions for the valence quarks. 
We calculate the two point function (2pt) II^p' and three point function (3pt) H^^', 



+,-5.(i,-.?o)^0|T [xN{l2,t2)J^{xuh)XN{xoM] |0), 



(2.1) 



where Xn is the nucleon interpolation field and the insertion is given by the point-split conserved 
current /^(x + /x/2) = (l/2)x s{x + + y^,)uI{x)s{x) - s{x){\ - Yfi)U ii{x)s{x + ^) . 

Electromagnetic form factors can be obtained using p = 0, p' = —^kinematics for the forward 
propagation (t2 S> fi ^ fo) l|5|]. In this work, we consider the backward propagation {t2 <C ?i <C to) 
as well, in order to increase statistics. The formulas for G|- and are summarized as 



<(rJoi) 



Tr 



^ pol 



Tr 



Tr[r±-n2pt(±^,?i; to) 



Tr 



r±-n2p'(o,fi; to) 



T^ijkqj 



pi . 



n2pt(o,?2; to) 



(2.2) 



(2.3) 



(±/)/2 X (1 ± 74)757^ and = Jmjj + q^. The upper 



where {i,j,k} / 4, r± = (1 ± 74)72 , r± 
(lower) sign corresponds to the forward (backward) propagation. Furthermore, we consider another 
kinematics of p = q, p' = 0, where the analogs of Eqs. ( 2^ ), ( 23 ) hold. We find that the results from 
the latter kinematics have similar size of statistical errors as those from the former, and the average 
of them yields better results. Hereafter, we present results from total average of two kinematics and 
forward/backward propagations, unless otherwise noted. In order to extract the matrix elements, 
we take the summation over the insertion time fi, symbolically given as 7?^ ^ = Lfi=fo+fj^^ ~ 



const . + f2 X G| ^ , where ^ are trivial kinematic factors appearing in Eq. (2.3) and t^ is chosen 
so that the error is minimal. We thus obtain G|. ^ as the linear slope of R'^ ^ against t2. 



The calculations of 3pt need the evaluation of DI. We use the stochastic method [15], with Z(4) 



noises in color, spin and space-time indices. We generate independent noises for different config- 
urations, in order to avoid possible auto-correlation. To reduce fluctuations, we use the charge 
conjugation and 75-hermiticity (CH), and parity symmetry [jll], |l^, |T7|]. We also perform unbiased 
subtractions [ [l8| ] to reduce the off-diagonal contaminations to the variance. For subtraction opera- 
tors, we employ those obtained from hopping parameter expansion (HPE) for the propagator 

where D denotes the Wilson-Dirac operator and C the clover 



2k 



M 



l+C ^ i+c' 



term. We subtract up to order (k"D)^ term, and the statistical error is reduced by a factor of 2. 

In the stochastic method, it is quite expensive to achieve a good signal to noise ratio (S/N) just 
by increasing Nnoise because S/N improves with y/N noise- In view of this, we use many nucleon point 



sources N^rc in the evaluation of the 2pt part for each configuration [17]. Since the calculations of 
the loop part and 2pt part are independent of each other, this is expected to be an efficient way. 
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Figure 1: R'^ (left) and R'j, (right) with K„d = 0. 13760, = 2- {27t /Laf, N,,rc = 64 (circles) and N,rc = 4 
(triangles, with offset for visibility), plotted against the nucleon sink time f?- The dashed hne is the linear fit 
where the slope corresponds to the form factor. 

In particular, for the Nnoise ^ Nsrc case, we observe that S/N improves almost ideally, by a factor 
of ^/N^■. We take A^„r = 64 for K^d = 0.13760 and A^,,.^ = 82 for K^d = 0.13800,0.13825, where 
locations of sources are taken so that they are separated in 4D-volume as much as possible. 



3. Results for strangeness form factors 

We calculate for the five smallest momentum-squared points, = n - {2n/La)^ {n = 0-4). 
Typical figures for R'^ ^ are shown in Fig. |l|. One can observe the significant S/N improvement 
by increasing Nsrc- Of particular interest is that, for all fc,,^/ simulations, is found to be 

negative with 2-3 a signals for low regions. In order to determine the magnetic moment, 
the dependence of G'l^{Q^) is studied. We employ the dipole form in the fit, Gl^{Q^) = 
^m(0)/(1 + 2^/A^)^, where reasonable agreement with lattice data is observed. For the electric 
form factor, we employ G|(e^) = g|r • Q^/il + Q^/A^ f, considering that G^(0) = from the 
vector current conservation. In the practical fit of G|-(2^), however, reliable extraction of the pole 
mass A is impossible because data are almost zero within error. Therefore, we assume that 

G|.(2^) has the same pole mass as G^(2^), and perform a one-parameter fit for JTl]]. 

Finally, we perform the chiral extrapolation for the fitted parameters. Since our quark masses 
are relatively heavy, we consider only the leading dependence on mk, which is obtained by heavy 
baryon chiral perturbation theory (HB;^PT). For the magnetic moment G^(0), we fit linearly in 
terms of mk [||, 19 1 . For the pole mass A, we take that the magnetic mean-square radius (r; 



2\ 



M 



12G^(0)/A^ behaves as 1/m^ For we use the electric radius {r^)E = —^Se which has 



an ln{mK/lJ.) behavior []T9|], and we take the scale /i = 1 GeV. The chiral extrapolated results are 
G^(0) = -0.017(25), Aa = 0.58(16), {rj)M = -7.4(71) x lO^^^fm^, {rj)E = -2.4(15) x lO^^fm^. 

Before quoting the final results, we consider the systematic uncertainties yet to be addressed. 
First, we analyze the ambiguity of dependence in form factors, by employing the monopole 
form [|ri|]. The results are found to be consistent with those from the dipole fit, and we take the 
difference as systematic uncertainties. Second, we study the uncertainties in chiral extrapolation 
by testing two alternative extrapolations. In the first one, we take into account the nucleon mass 
dependence on the quark mass, using the lattice nucleon mass. From the physical viewpoint, this 
corresponds to measuring the magnetic moment not in units of lattice magneton but physical mag- 
neton [pO|]. In the second alternative, we use the linear fit in terms of m^, observing the results have 
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Figure 2: The chiral extrapolated results for G^(2^) (left) and Gg{Q^) (right) plotted with solid lines. 
Shaded regions represent the error-band with statistical and systematic error added in quadrature. Shown 
together ai-e the lattice data (and g^-extrapolated G^(0)) for = 0.13760 (circles), 0.13800 (triangles), 
0.13825 (squares) with offset for visibility. 



weak quark mass dependence. In either of alternative analyses, we find that the results are consis- 
tent with previous ones. While a further clarification with physically light quark mass simulation 
and a check on convergence of HB;^PT [ 21 ] is desirable, we use the dependence of results on differ- 
ent extrapolations as systematic uncertainties. Third, we examine the contamination from excited 
states. Because our spectroscopy study indicates that the mass of Roper resonance is massive com- 
pared to the ^ii state on the current lattice [22], the dominant contaminations are (transition) form 
factors associated with Si i . On this point, we find that such contaminations can be eliminated the- 
oretically, making the appropriate substitutions for Ff in Eq. (22) and {Tf, F^} in Eq. (23) [jlll]. 
It is found that the results from this formulation are basically the same as before, so we conclude 
that the contamination regarding the i state is negligible. 

As remaining sources of systematic error, one might worry that the finite volume artifact could 
be substantial considering that the spacial size of the lattice is about (2fm)^. However, we recall 
that Sachs radii are found to be quite small, \{r^)E,M\ O.lfm^, which indicates a small finite 
volume artifact. For the discretization error, we conclude that finite {qa) discretization error is 
negligible, since the lattice nucleon energy is found to be consistent with the dispersion relation. 
As another discretization error, we note that {mK) is found to have 6 (8) % error for the current 
configurations [14, 23]. Considering the dependence of ^ on these masses, we estimate that the 
discretization errors amount to ^ 10%, and are much smaller than the statistical errors. Of course, 
more quantitative investigations are desirable, and such work is in progress. 

To summarize the results of form factors, we obtain G^(0) = —0.017(25) (07), where the 
first error is statistical and the second is systematic from uncertainties of the extrapolation and 
chiral extrapolation. We also obtain Aa = 0.58(16)(19) for dipole mass or Ka = 0.34(17)(11) 
for monopole mass, and g| = 0.027(I6)(08). These lead to G'^iQ^) = -0.015(23), ^^(2^) = 
0.0022(19) at = O.lGeV^ where error is obtained by quadrature from statistical and systematic 
errors. We also obtained, e.g., G'^iQ^) = -0.014(21), G^Q^) = 0.0041(38) at = 0.22GeV2. 
Note that these are consistent with the world averag ed data at = O.lGeV^ [§ g |] and the 
recent measurement at Mainz [|4|], 0^(2^) = -0.14(11)(11), G^Q^) = 0.050(38)(19) at = 
0.22GeV^, with an order of magnitude smaller error. In Fig. ^, we plot our results for G^(2^), 
G|.(Q^), where the shaded regions correspond to the square-summed error. 
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Figure 3: The ratio of 3pt to 2pt with K^d = 0.13760, = {2n/La)^, A(,„- = 64 (circles) and N,rc = 4 
(triangles, with offset for visibility), plotted against the nucleon sink time t2- The dashed line is the linear fit 
where the slope corresponds to the second moment. 



4. Results for the strangeness second moment {x^)s- 



In the calculation of the strangeness second moment, {x^)s-s, we use the following three-index 
operator as the insertion operator in 3pt, 



'4(7 



1 

'3 L' 



D i D iS + sYi D 4 D jS + sji D , D 45 



(4.1) 



Taking the kinematics of p' 
moment (ITi 



p,q = 0, the following ratio of 3pt to 2pt corresponds to the second 



Tr 



r±-ng(±^,f2; OA; ±a?o) 



Tr[r±-n2pt(±^,f2; to) 



(4.2) 



where the upper (lower) sign corresponds to the forward (backward) propagation as before. 

In the evaluation of the DI, we can apply basically the same technique which is used in the 
calculation of the strange form factors. Because of the difference of the structure of the insertion 
operator, we perform the unbiased subtraction from HPE up to order {kD)^ term. 

In Fig. |3| , we plot the ratio of 3pt to 2pt in terms of the nucleon sink time, t2, where the linear 
slope corresponds to the signal of {x^)s-s- One can clearly observe that S/N improves significantly 
by increasing Ns,-c, as was observed in the form factor study. Yet, the signal of {x^)s-s in this figure 
is still consistent with zero. In order to obtain the final result, detailed analyses are in progress. 



5. Summary 

We have studied the strangeness electromagnetic form factors of the nucleon from the Nf = 
2+1 clover fermion lattice QCD calculation. It has been found that calculating many nucleon 
sources is essential to achieve a good S/N in the evaluation of DI. We have obtained the form 
factors which are consistent with experimental values, and have an order of magnitude smaller 
error. The importance of the strangeness second moment, {x^)s-s, has been emphasized, and a 
preliminary result has been reported. 
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